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Abstract 

Self-organization of proteins in space and time is of crucial importance for the functioning 
of cellular processes. Often, this organization takes place in the presence of strong random 
fluctuations due to the small number of molecules involved. We report on stochastic switching 
of the Min-protein distributions between the two cell-halves in short Escherichia coli cells. 
A computational model provides strong evidence that the macroscopic switching is rooted in 
microscopic noise on the molecular scale. In longer bacteria, the switching turns into regular 
oscillations that are required for positioning of the division plane. As the pattern becomes 
more regular, cell-to-cell variability also lessens, indicating cell length-dependent regulation of 
Min-protein activity. 
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I. INTRODUCTION 



Subcellular structures are often formed by a small number of proteins. This is notably the case 
in prokaryotes, as exemplified by the Soiproteins in Bacillus subtilis that switch randomly between 
the two distal parts of the nucleoid [l[ Q], by the Min proteins in Escherichia coli jij] that switch 
periodically between the two halves of the bacterium, as well as by the rings and helices formed by 
FtsZ and MreB in many bacteria [4|. The concentrations of these proteins are in the micromolar 
range, hence, the above structures, that may extend over several micrometers, are formed by a 
few hundred molecules. Such small numbers can imply large random fluctuations in space and 
time 

However, there is currently a lack of quantitative experimental studies of random fluctuations 
in spatially extended protein structures. In particular it remains largely unexplored how fluctua- 
tions on cellular scales originate from the interplay of molecular noise, i.e., the stochastic nature 
of the protein kinetics, and external noise, i.e., variations in the environment or in the number 
of proteins present in the cell 4-t| Q|- In the context of gene expression, random fluctuations 
contribute essentially to the dynamics [l(J. In particular, they directly lead to differences in ge- 



Ill ] . i.e., intercellular fluctuations, and can play an important role in cell 



netically identical cells 

fate decisions [ill. \l2\. In this context . cell size can act as a control parameter to regulate either 
the amplitude of fluctuations [13|] or the distribution of phenotypes they generate 



HQ- 



In this work, we report on intra- and intercellular fluctuations in the spatiotemporal patterns 
formed by the Min proteins in E. coli. These proteins select the cell center as the site of division [jj] . 
Their distribution has been found to change periodically with time, such that most of the proteins 
reside for about 40s in one cell half and subsequently for the same time in the opposite cell half fis| . 
As one of the Min proteins, MinC, inhibits formation of the Z-ring, division is in this way suppressed 
at the cell poles. Computational analysis of the Min-protein dynamics indicate that the observed 
pattern results from self-organization of the ATPase MinD and MinE [16]. This idea is supported 
by experiments in vitro in which planar and spiral waves of MinD and MinE emerged spontaneously 
in the presence of ATP on a supported lipid bilayer 

We found, that in cells below a critical length of 2.7//m the Min distribution switches stochas- 
tically between the two cell halves. In this phase intercellular fluctuations are pronounced. In 
cells longer than the critical length, we observed regular oscillations and intercellular fluctuations 
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were significantly reduced. Computational analysis shows that stochastic switching can result 
from molecular noise. Furthermore, it indicates that, in the stochastic state, the rate of ATP 
consumption per MinD molecule is only half of that in the oscillatory state. 

II. RESULTS 

We observed the distribution of MinD for 40 minutes in each of 209 cells, see Materials and 
Methods, during which individual bacteria grew roughly 0.5/im. In cells shorter than 2.5/xm, instead 
of oscillating regularly, MinD typically shifted stochastically from one cell half to the other, see 
Fig. [Th.. The residence times of MinD in one cell half varied widely in these cells, whereas complete 
switching from one cell half to the other occurred in an interval of less than 15s. The extent of the 
region covered by MinD in one cell half did not change notably between two exchange events, see 

Fig. dg. 

For cell lengths between 2.5/im and 3^m, the Min pattern typically changed from stochastic 
switching to regular oscillations with a period of about 80s, see Fig. [lb and Fig. SI in the Supp. Mat. 
The precise lengths at which this transition occurred differed between cells. In the oscillatory 
regime, between two switching events, the region covered by MinD first grew monotonically from 
the cell end and then shrank monotonically towards the same end, see Fig. [Tg. Analogously to the 
stochastic exchange regime, the transition of MinD from one half to the other was fast compared 
to the oscillation period. In a small number of cells (N=5), we observed a transition back from 
regular oscillations to stochastic switching. 

Cells longer than 3.5/im invariably displayed regular oscillations, see Fig. [lb and Fig. SI in the 
Supp. Mat. The oscillation period typically decreased slightly with increasing cell length. For the 
cell shown in Fig. [Tb, the initial period is approximately 87s, while it is approximately 70s at the 
end. Similar behavior can be observed for MinE, see Supplementary Material. 

About 5% of the cells divided during the observation time, see Fig. [TH, e. In all these cases, 
MinD oscillated regularly prior to division. Consistent with our findings in non-dividing cells, the 
pattern displayed by the Min proteins in the daughter cells immediately after division correlated 
with their length: daughter cells shorter than 2.5/um typically displayed stochastic switching, while 
the Min proteins mostly oscillated in daughter cells longer than 3/im. In some cases, however, 
the two daughter cells showed different MinD patterns, in spite of having equal lengths: while in 
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one daughter we observed oscillations, the other displayed stochastic shifts. In these cases, judged 
from the fluorescence intensities, the distribution of MinD between the two daughter cells was 
significantly uneven 31]. 



A. Distribution of residence times 

The transition from stochastic switching to regular oscillations is accompanied by a qualitative 
change in the distribution of residence times r cumulated over cells of a given length, see Fig. [2J 
For cells longer than 3//m, this distribution shows a pronounced peak at ~ 35s, which essentially 
coincides with the mean value, and a small tail including rare events of residence times of up to 
200s. These rare events correspond to instances when one or more oscillations were missed and 
thus the residence time extends to values significantly larger than the mean, see Fig. [9] in the 
Supplementary Material. The left, main part of the distribution (r < 60s) can be fitted by a 
log-normal distribution with a geometric mean of 37.9s and a geometric standard deviation of 1.2s. 
The tail (r > 60s) can be fitted to an algebraic decay with decay exponent a = 3.6 ± 0.6. 

For cells smaller than 2.5/im, the distribution of residence times is qualitatively different: it is 
essentially algebraic, oc r _a , with a decay exponent a = 2.1 ± 0.2. Consequently, the mean-value 
is just at the edge of being well-defined, while the variance diverges, implying enormous variations. 
This is in contrast to a usual random telegraph process in which a system switches stochastically 
between two states at given constant rates. In this case, the distribution of residence times in 
one state decays exponentially. We, therefore, also fitted an exponential decay to the tail of the 
distribution of residence times. However, the logarithm of the likelihood ratio of an algebraic to 
an exponential decay was 35 and thus strongly indicative of an algebraic decay fis| . 

In Fig. we show the functional dependence on the cell length of the mean residence time 
(r) cumulated over all events and all cells in a given length interval. Cell lengths were binned into 
intervals of 0.2/im length and residence times were assigned to the length of the cell at the beginning 
of the respective residence time interval. Note, that since within 500s, which is exceptionally 
large for a residence time, the bacteria grew at most 0.2/xm in length, our results do not depend 
significantly on the exact assignment rule. Strikingly, the variation of (r) is well-described by 
two exponential decreases: For cell lengths smaller than 2.7//m that are typically in the stochastic 
regime, (r) cx exp(— L/X) with A = 0.56^m, while for longer, typically oscillating cells A = 2/iin. 
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To characterize fluctuations around (r). we computed the effective standard deviation a corre- 



sponding to our finite number of samples [32j. The effective standard deviation <?, too, decreases 
exponentially with increasing cell length. Note that, for cells smaller than approximately 2.7/iin 
it is larger than (t), while for longer cells it is smaller. For cells larger than 3.5/im, the standard 
deviation is below 10s indicating that the oscillation period varies remarkably little within a single 
cell and within an ensemble of cells. 



B. Cell-cell variability 

To assess the contribution of cell-to-cell variability to the fluctuations, we compared the mean 
residence time fi c of individual cells with that of other cells sharing the same cell length. In 
Fig. [3b, we represent the normalized standard deviation of the distribution of mean residence 
times, o"i c = \J (t?,) — (ilc) 2 / (t\c)- It decreases with increasing cell length, first mildly then sharply 
for cell lengths larger than 3pm, indicating a clear reduction of cell-to-cell variability in longer 
cells. This suggests a control mechanism that adjusts Min parameters towards optimal reference 
values when the cells grow. 



C. Mathematical modeling 



Several attempts at mathematically modeling of the Min system have been successful in re- 
producing essential aspects of the Min dynamics |lg]. Based on different underlying microscopic 
pictures, these models uniformly explain the Min oscillations as an emergent property of inter- 
acting MinD and MinE in presence of a membrane and ATP. In particular, self-enhanced binding 
of MinD and/or MinE to the membrane, or aggregation of membrane-bound proteins can trigger 
an instability towards oscillations. Biochemically, MinD and MinE are not yet characterized well 
enough to exclude one possibility or the other. The validity of a model can alternatively be tested 
by verifying predictions of the macroscopic behavior. For example, the aggregation models in 



Refs. 



19, 



2Cj] predict the existence of stationary heterogeneous protein distri 
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nutions. 



201 ] and checked if 



We tested several stochastic versions of models presented in the past 
they showed stochastic switching of Min proteins in certain parameter regimes, see Supp. Mat. 
Out of these, only the model suggested in [20J produced stochastic switching similar to what is 
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observed in the cell 331. Since stochastic simulations in three dimensions yield essentially the same 

Q Q 

results as in one dimension, cf. Refs. [7[ and [9(, we restricted our attention to simulations in one 
dimension. 

What is at the origin of the transition between stochastic switching and oscillations of MinD 
and the corresponding change in the residence-time distribution? The experimental data presented 
above show that the spatiotemporal structure of the MinD pattern correlates with cell length. As 
a naive guess one might think that the cell length directly controls the pattern. However, in our 
experiments, we occasionally observed that two daughter cells of about the same length showed 
different MinD patterns right after division of the mother cell. In addition, in some instances, cells 
shorter than 2.5/im also showed oscillations, while in some cells longer than 3fim MinD switched 
stochastically between the two cell halves. This suggests that, in addition to the cell length, other 
factors control the transition from stochastic switching to regular oscillations. 

In order to better understand which parameters of the system are able to trigger stochasticity in 
Min switching, we used a one-dimensional particle-based stochastic version of the one model which 
we identified as being able to produce stochastic switching, see Fig. [Hand Material and Methods. 
We find that, increasing the cell length is not sufficient to generate a transition from stochastic 
switching to oscillatory Min dynamics if, simultaneously, protein numbers increase proportionally 
to the cell length. Changing other parameters of the Min dynamics, however, can trigger such 
a transition. As an example we now discuss alterations in the rate loe of MinE binding to the 
membrane, see Fig. [5j Let us stress that in the following we keep the total protein concentrations 
constant, that is, protein numbers vary proportionally to the cell length. For ue = 0.1s -1 and a 
cell length of 2/im, our simulations show stochastic switching, while for loe = 0.35s -1 and a cell 
length of 3/xm, we find oscillations 



34 ] , see Fig. |5k. The model is thus capable of reproducing the 
observed transition by changing MinE activity. It is important to note, that this transition is not 
a consequence of approaching the deterministic limit, but is inherent to the Min-protein dynamics. 



Indeed, for the respective binding rates, the deterministic model is either bistab 
one of two mirror-symmetric stationary states or it generates oscillations 



3, 



e and settling into 
2o| . We conclude 



that in short cells the macroscopic switching behavior is rooted in intrinsic fluctuations of the 
molecular processes due to the relatively small number of proteins involved. The transition from 



stochastic switching to oscil 
which in general is possible 



ations, however, is not due to the discreteness of molecule numbers, 
2J. 



Quantitative analysis of simulations corresponding to a single cell in the stochastic switching 
state, presented in Fig(5k, yield an exponential decay of the residence times as expected for a ran- 
dom telegraph process. This is however not in contradiction with the experimental data, for which 
ensemble averaging over different cells yields algebraic decay, Fig. [2h.. Indeed, experiments revealed 
cell-to-cell variability, particularly strong amid short cells exhibiting stochastic switching, Fig. [3b. 
As we show below, introducing such variability in the model resolves the apparent contradiction. 
It is conceivable that in the experiments the cell-to-cell variations in the Min dynamics for a given 
cell length result mainly from differences in the numbers of MinD and MinE. Indeed, we found in 
~ 65% of the observed divisions that the numbers in the daughter cells, as judged from fluorescence 
intensities, differed by more than 10%, see Supplementary Material. 

In order to test the effects of varying protein numbers in the model, we performed several 
simulation runs with MinD protein numbers drawn at random from a Gaussian distribution with 

n 



a standard deviation of 10% of the mean [ill]. The ratio of MinD to MinE was fixed to 8/3 [2] and 
all other parameters kept constant. Lumping together the observed residence times of 70 runs with 
different protein numbers, we now observe an algebraic decay with exponent 2.03 ±0.2, see Fig. [5b, 
which is very close to the experimental value of 2.1, see Fig. [2k.. As in the experiments, the mean 
is barely defined and the standard deviation diverges. This result constitutes a second instance 
of the general finding that power law distributions can result from Gaussian variations of s yst em 



25] 



parameters as discussed by Tu and Grinstein in the context of the bacterial flagellar motor 

Analogously, we calculated the distribution of residence times in the oscillatory regime, see 
Fig. 03. It is remarkably similar to that found experimentally, see Fig. [2b: a large hump, well 
fitted to a log normal distribution with geometric mean 31s and geometric standard deviation 1.2s, 
plus a small tail at large residence times which can be fitted by an algebraic decay with exponent 
a = 4.6 ±0.4, similar to the experimental value a = 3.6. Moreover, an inspection of the rare events 
at the origin of this algebraic tail showed that they are instances of "missed" oscillations, just like 
in the experimental kymographs, see Fig. [9] in the Supplementary Material. We conclude that the 
oscillation period is a robust property of the Min System with respect to fluctuations in protein 
numbers. 

We then investigated the transition from the switching to the oscillatory state in more detail. 
Fig. [5H presents the mean value and the standard deviation as a function of system length. The 
rate loe was varied with the system length according to a sigmoidal dependence, see Fig. [Sfe inset 



and Supplementary Material. Similar to the experimental results, we find a decreasing mean 
and standard deviation with increasing length and the transition from stochastic switching to 
oscillations is accompanied by a drop of the standard deviation below the mean value. The simple 
particle model of the Min dynamics thus semi-quantitatively captures the effects of length changes 
and of fluctuations in the system. 

Finally, we used the computational model to infer the ATP consumption rate as a function of cell 
length. Interestingly, the average ATP consumption rate per unit length, which is proportional to 
the ATP consumption rate per MinD, increases with cell length until the transition from stochastic 
switching to regular oscillations and then remains roughly constant, see Fig. [5fe. Our calculations 
thus let us hypothesize that is energetically advantageous for a short cell to keep the Min system 
in the stochastic switching regime. 



III. DISCUSSION 

In summary, we have shown that in short E. coli cells the distribution of MinD proteins stochas- 
tically switches between two mirror-symmetric states, while it oscillates regularly in longer cells. 
The Min system provides thus an intriguing example of a spatiotemporal pattern under physiolog- 
ical conditions that combines regularity with stochastic elements. 

Stochastic switching between two states is also known for the Soj/SpoOJ proteins which in 
B. subtilis relocates irregularly between the two sides of the bacterial nucleoid [l|, |2|. A computa- 
tional model suggests that this phenomenon, too, is rooted in molecular fluctuations [6(. Similarities 
between MinD and Soj had already been pointed out due to resemblance in structure, polymeriza- 
tion and ATPase activity j^. However, notable differences exist as, for example, SpoOJ seems to 
be always bound to the nucleoid unlike MinE that detaches from the membrane. Furthermore, no 
oscillations have been found for Soj/SpoOJ. 

Based on our theoretical analysis, we can speculate about a possible regulatory mechanism con- 
trolling the dynamics of the Min proteins in E. coli. We interpret the findings of our computations 
as indications of changes in the activity of the Min proteins as the cells grow. Specifically, the 
ability of MinE to bind to membrane bound MinD might be reduced in early phases of the cell 
cycle. According to this view, the Min system is kept in a stand-by mode with reduced ATP 
consumption in short cells and gets fully activated only as a cell approaches division, when the 
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Min system needs to be functional. Why does stochastic switching of the Min proteins not lead to 
minicelling? The central Z-ring needs about 5 to 10 minutes to mature and polar Z-rings are less 
stable than the central one. Consequently, polar rings formed in an early phase of the cell cycle 
might be disassembled by MinC after oscillations have set in, which occurs significantly before 
septation starts. 

Obviously, the transition might be also due to other evolutionary constraints or be just a side- 
effect of the Min-protein dynamics. Independently of the transition's physiological reason, we 
propose that changing system parameters such that the dynamic behavior is qualitatively modified 
through a bifurcation - in the present case from bistable to oscillatory - presents an interesting 
generic mechanism to tame detrimental fluctuations as subcellular processes become vital. 

Materials 

A. Data acquisition 

We used cells of the E. coli strain JS964 containing the plasmid pAM238 encoding for MinE 
and GFP-MinD under the control of the lac-Promoter (2^]. Bacteria were grown overnight in a 
3ml LB medium at 37° C. Cells were induced with Isopropyl-/3-D-thiogalactopyranosid (IPTG) at 
a concentration of 200/uM and incubated for 4 hours. During measurements, cells were in the 
exponential growth phase. The samples were kept at a temperature of 30° C using a Bachhoffer 
chamber. To keep bacteria from moving under the cover slip, we put them on an agar pad (1% 
agar solution in LB medium with a reduced yeast extract fraction, 10%, in order to lower back- 
ground fluorescence). The fluorescence recordings were taken with an Olympus FV 1000 confocal 
microscope, at an excitation wavelength of 488nm from a helium laser at low power. We used 
an Olympus UPLSAPO 60x, NA 1.35 oil immersion objective and recorded a frame every 3s. A 
measurement lasted 40min. During this period, the focus was manually readjusted at irregular 
intervals. We were not able to determine the protein numbers in individual cells. 

In total, we extracted data from 209 MinD-fluorescent cells obtained from 5 different measure- 
ments. We extracted data only for cells which, at the beginning of the measurement, were smaller 
than 3/iin. Cell lengths were determined from differential interference contrast (DIC) images with 
an accuracy of ±150nm at the beginning and the end of a measurement. The cell length in-between 
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was determined by linear interpolation. Some of the cells in the field of view divided during the 
measurement time. If division occurred after more than 20min of measurement, they were in- 
cluded in the data analysis. In these cases, fluorescence recordings were used until cell constriction 
terminated. 

To observe MinE distribution in cells, we used the strain WM1079 expressing MinD and MinE- 
GFP on a plasmid under the control of the pBAD-promoter [l|. Bacteria were induced with 
Arabinose but otherwise grown in the same way as described above. 



B. Data analysis 

For a quantitative analysis, we mapped the time-lapse data of MinD-fluorescence from E. coli 
cells onto the states of a two-state process. We transformed the fluorescence data into a real- valued 
time-dependent function by subtracting the integrated fluorescence intensities from the left and 
the right cell half, see Fig. Then, the moving average over four time points was taken in order 
to reduce noise. The resultant function / is positive when the fluorescence maximum is in one cell 
half and negative in the opposite case. The residence time r of MinD in one cell half is defined 
as the interval between two consecutive zeros of /. To each residence time, we assigned the cell 
length at the start of the respective residence period. 



C. Computational model of the Min dynamics 



For the computational ana" 
of the model suggested in 



lysis of the Min-protein dynamics, we use a particle-based version 
20 ]. It describes the formation of Min oscillations on the basis of 
an aggregation current of bound MinD that is generated by mutual attraction of the proteins. 
Furthermore, it accounts for the exchange of MinD and MinE between the membrane and the 
cytoplasm, where MinE only binds to membrane-bound MinD and where MinD detaches from 
the membrane only in the presence of MinE. Measurements of the cytoplasmic diffusion constants 
of MinD and MinE yielded values larger than 10/xm 2 /s We thus consider the limit of large 
cytoplasmic diffusion which effectively leads to homogeneous cytoplasmic concentrations. The 
diffusion constant of membrane-bound proteins is about two orders of magnitude smaller than for 
cytosolic diffusion 29]. This mobility suffices to generate a sufficiently strong aggregation current 
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generating an instability of the homogenous protein distribution 20] . 

A corresponding stochastic particle-based description is defined on a one-dimensional lattice 
with N sitesrepresenting the long axis of the cell. The lattice spacing If, is chosen such that it 
is much larger than the protein size and much smaller than the characteristic length of the Min 
pattern. Each site can contain at most n max proteins, a number which can be understood as the 
circumference of the cell. We assume diffusional mixing such that proteins are indistinguishable on 
a site. For site j, the probability of attachment of cytoplasmic MinD and MinE during a sufficiently 
small time step At is given by 



P D ^ d = At u D (1 - ' (1) 

\ iv / ""max 

P E ^ = AtuJ^)^L , (2) 

respectively. Here, ujd,oje are the corresponding attachment rates and Njj and Ne are, respec- 
tively, the numbers of cytoplasmic MinD and MinE. The numbers of membrane-bound MinD and 
MinDE complexes on site j are n d j and n de j. The detachment probability is 

Pde^E+D = At LO de n de j , (3) 

where oj de is the detachment rate. The exchange of particles between sites is governed by 

D Dd^t n d,j±i + n de j±i 

Pj->j±l = — 72— n d ,j{l — J —)Ij^ j±1 , (4) 

"max 

where 

f 1 if AEj < 

^•±l = \ AF ( 5 ) 

[exp(-0) if AE,>0 

with AEj = Vj±i — Vj. The potential V describes the interaction strength between Min-proteins 
on the membrane. We assume a square hole potential 



(2R d + 1) P d + (2R de + 1) ,P o ^ 



(6) 



j = ~Rd j = -Rde 

Here, the integers R d and R de relate to the ranges r d of the MinD-MinD interaction and r de of the 
MinD-MinDE interaction through R d ~ Td/h and Rde — r de/h- The parameters g d and g de tune 
the interaction strength. The diffusion constant of membrane-bound MinD is D d . 

In Fig. [5]i, the attachment rate is increased jointly with the cell length in the simulations. 

11 



We chose a sigmoidal increase of wg in dependence of the cell length £ c according to the Goldbeter- 
Koshland function, which gives the mole fraction of modified proteins that are under control of a 
modifying enzyme j^] 

G(4) = -»x2^ (7) 

vJ + £ C (K - 1) + y/(vJ + i c {K - l)) 2 - 4(u - l c )l c K 

We chose parameters v = 2.47, J = 1.116 and K = 0.099. The factor ujg determines the saturation 
value and was chosen as 0.4s -1 . The precise functional form of oje{&c)i however, is not of importance 
as long as it is sigmoidal. 
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FIG. 1: Distribution of MinD-GFP along the bacterial long axis as a function of time. Time 
increases to the right, the long axes are oriented vertically. Initial cell lengths are given on the 
left, final cell lengths on the right, a) Stochastic exchange of MinD between the two cell halves, 
b) Stochastic exchange turned into regular oscillations at a length of 2.8/Uin (green arrow), c) 
Regular Min oscillations with a period of 73s. d, e) Min dynamics in dividing cells. Daughter cells 
show either stochastic switching (d) or regular oscillations (e). Divisions occurred, respectively, at 
lengths of 4/im and 5/im of the mother cells (red arrows), f) Nomarsky and fluorescence images of 
the cells used in (a)-(e) at the beginning of the measurements. In (f5), lines indicate the bacterial 
long axis and width, g) Extracts of the kymographs shown in a) and c), respectively. 
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FIG. 2: Ensemble distributions of residence times, a) Residence times for initial cell lengths 
between 2 and 2.2\xm. (total: 358) and b) for initial lengths greater than 3/xm (total: 2157). Times 
are respectively binned into 20s (a) and into 3s intervals (b). Solid lines result from fits of the 
distribution tail (r < 60s), (a) Fitting of an algebraic tail yields an exponent —2.1 ± 0.2 (red line), 
(b) the main part of the distribution is well fitted to a log-normal distribution with geometric mean 
37.5s and geometric standard deviation 1.2s (blue line); the tail of rare events of duration longer 
than 60s is reasonably well fitted to an algebraic law with exponent —3.6 ± 0.6 (red line). 
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FIG. 3: a) Mean (r) and effective standard deviation a of ensemble distribution of residence times 
as a function of cell length L. Cell lengths have been binned into 0.2^m intervals. Straight lines 
are exponential fits to the (r) data, r = roexp(— L/X). For cell lengths below 2.7/im we find 
A = 0.5/im, for larger lengths A = 2/xm. b) Relative standard deviation a lc = J (f? c ) — (r ic ) 2 / (r ic ) 
of the distribution of individual cell mean residence times 

lie cl function of cell length. 
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FIG. 4: Stochastic model of Min-protein dynamics. The cell membrane is represented by a one- 
dimensional array of bins that can take up MinD-ATP (green) and MinE (red). Cytosolic MinD- 
ATP binds at rate ujd, cytosolic MinE binds to membrane-bound MinD-ATP at rate uje, MinD- 
ATP-MinE complexes detach at rate uJde- The distribution of proteins in the cytosol is assumed 
to be homogenous [29I ] and the exchange of MinD-ADP to MinD-ATP infinitely fast. On the 
membrane MinD-ATP diffuses and is attracted towards other MinD-ATP. For details see Materials 
and Methods. 
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FIG. 5: Theoretical analysis of the Min-protein dynamics, a) Kymographs obtained from sim- 
ulations with ue = 0.1s" 1 and length 2/xm (top) and wg = 0.35s -1 and length 3/xm (bottom), 
other parameters as below. Scale bar indicates 5min. b,c) Ensemble distribution of residence 
times for the case of stochastic switching (wg = 0.1s -1 , cell length: 2/im.) and oscillatory dynamics 
(oje = 0.35s -1 , cell length:3/im). Each histogram was obtained from 70 runs where the number 
of MinD was drawn from Gaussian distributions with mean 1440 (b) and 2160 (c), respectively 
and 10% standard deviation. The ratio of MinD to MinE proteins was fixed at 8/3. The red lines 
represent algebraic fits r~ a for residence times r > 60s. For (b) a = 2.03±0.2, for (c) a = 4.6±0.4. 
In (c), the blue line indicates a log-normal distribution with geometric mean 31s and geometric 
standard deviation 1.2s. d) Dependence of the mean residence time (r) (green circles) and effective 
standard deviation a (black squares) on the cell length obtained from ensemble simulations as in 
(b) and (c). e) The rate of ATP hydrolysis per unit length as a function of cell length. Inset: 
dependence of the MinE attachment rate uje on cell length. Other parameters are oje = 0.04s -1 , 
UJde = 0.04s -1 , D d = 0.06/um 2 /s, r d = 1.2/im, r de = 0.1/iin, g d = 35k B T, g de = -20k B T, n max = 43 
and bin length Z& = 33nm. 
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FIG. 6: Kymograph of GFP-MinD fluorescence and corresponding intensity curve / obtained by 
subtracting the fluorescence intensity of one cell half from the intensity of the other cell half. 
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IV. SUPPLEMENTARY MATERIAL 



A. Stochasticity of Min Switching Decreases in Individual Cells During Cell Elon- 
gation / Aging 

We analyzed the stochasticity of MinD switching in individual cells for time intervals in which 
the cell grew less than Al. For the associated Min dynamics in such an interval, we calculated the 
ratio r lc of the standard deviation a lc of the residence times and the corresponding mean fi c , 

Tic — "lc/lc- 

This quantity measures the stochasticity of the Min dynamics in an individual cell - the higher it 
is, the more stochastic is the switching behavior. In Fig. we show a histogram of the number 
of cells which switch stochastically/regularly in a length interval of Al = 0.5^m according to the 
stochasticity measure n c . 

We averaged the quantity r lc over ensembles of cells in a cell length interval Al = 0.2/im. 
In Fig. [HJ we present the average ratio (rj c ) as a function of the cell length, (yellow triangles), 
where (...) denotes the average over all measured cells in the respective length interval. To get a 
meaningful estimate of r- lc , only cells showing at least 5 switching events within a length interval 
of 0.2//m were considered. Very long residence times are thus neglected, which, together with the 
small number of events entering n c , tends to systematically underestimate the standard deviation 
of the "real" distribution. The value of r initially drops monotonically and stays constant for cell 
lengths larger than 3.5//m. The ratio of the mean residence time and the corresponding standard 
deviation decreases thus in individual cells. 

For comparison, we took the residence times of the same cells and pooled them directly according 
to the appropriate cell length interval. This means that the averaging over cell populations is 
done first, before any statistical analysis. From the obtained sets of residence time distributions, 
we calculated again standard deviation a and mean (r) and their ratio r = <t/(t). Again, this 
quantity r drops for increasing cell length, Fig. [5] (green triangles), but takes on average higher 
values than (r; c ) (yellow triangles). This is because the standard deviation results in this case also 
from cell-cell variability. 
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B. Regularly Oscillating Cells, Protein Translocation Events Are Occasionally 
Missed 

Rarely, in regularly oscillating cells, a switching event is missed and the fluorescence maximum 
stays longer in one cell half. An example of such an event is shown in the top panel of Fig. [9j This 
phenomenon can also be observed in the simulated Min oscillations as shown in the bottom panel 
of Fig. [9l These events give rise to a small tail in the residence time distribution of oscillating cells. 

C. MinE Ring Stochastically Switches in Short Cells 

During our measurements of MinD fluorescence in cells, we could not assess the distribution 
of MinE simultaneously. Instead, we observed the distribution of fluorescently labeled MinE in 
a different Escherichia coli strain, WM1079 (SI). This strain produces on average larger cells 
and does not show a clear transition from stochastic Min switching to Min oscillations during cell 
growth. We did however observe stochastic switching of MinE in individual cells, Fig. [TO) This fact 
together with known results of simultaneously recorded distributions of MinD and MinE in vivo 
(S2) and in vitro (S3) suggests strongly that MinE is in the same way distributed in stochastically 
switching cells as in normally oscillating cells: A MinE ring at the rim of the MinD tube and a 
shallow MinE layer on the remainder of the MinD tube. 

D. MinD is Partitioned Unevenly Between the Daughter Cells During Cell Divi- 
sion 

We measured the total fluorescence intensity from daughter cells after cell division (the estimated 
error is about 10%). We then calculated the associated deviation of the daughter fluorescence from 
the equipartition value, i.e., the total fluorescence intensity which would be expected if MinD was 
partitioned equally to both daughter cells. Data from 37 cell divisions give a distribution of relative 
deviations from the equipartition value with a standard deviation of 21%, Fig. [TTJ The distribution 
peaks in the center meaning that a zero deviation is most likely. This is in contrast to the Min 
distribution suggested by the theoretical results presented by Tostevin and Howard (S4). There, 
the authors probed the partitioning of Min proteins using a one-dimensional stochastic model 
including MinD polymerization and the formation of MinD/MinE complexes on the membrane. 
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They simulated the division process assuming a gradually decreased diffusion constant in the cell 
middle due to cell constriction. As a result, they predicted a distribution of MinD fractions in the 
daughter cells which peaked at 20% deviation from the equipartition value. 

E. Increasing MinE Attachment Rate can Generate a Transition from Stochastic 
Switching to Regular Oscillations in Cells of Fixed Length 

We tested how the Min-protein dynamics changes within our stochastic one-dimensional model 
with an increasing MinE attachment rate loe in cells of a fixed length of 2.5/xm. We chose a 
mean MinD protein number of 1800 and a fixed MinD/MinE ratio of 8/3. For each value of 
we we performed 70 simulations with simulation time 160min each. To mimic the variation of 
actual protein numbers in real cells, we drew the protein numbers of MinD and MinE from a 
Gaussian distribution with 10% standard deviation. The remaining parameters of the simulation 
were wo = 0.04s -1 , = 0.04s -1 , = 0.06/xm 2 /s, = 1.2/xm, r^ e = O.l^m, = 35&bT, 
gde = —20ksT, n max = 43 and bin length lb = 33nm. As a result, we find that the stochasticity 
of the switching is reduced for increasing uje coupled with a drop in mean residence time, Fig. 1121 
For uje > 0.3s -1 we observe regular oscillations. 

F. Stochastic Model Switching Can be Triggered by Other Parameters Than uje- 

As mentioned in the main text and in the previous paragraph, our model predicts a transition 
from stochastic Min switching to regular oscillations by an increase of the binding rate of MinE 
to the cytoplasmic membrane. In addition, other model parameters can trigger such a transition 
in the dynamics (S5). It can be induced by i) an increase of the MinE/MinD ratio in the cell, 
by ii) a decrease in the density of binding sites on the membrane n max , by iii) a joint increase of 
MinE and MinD concentrations, and by iv) a change in MinD binding or MinD/MinE detachment 
from the membrane (data not shown). A pure increase of the cell length keeping the remaining 
system parameters constant (including protein concentrations, not protein numbers) results in a 
trend towards more stochastic switching for longer cells. This is the contrary effect to what we 
observe experimentally. This effect can be (over)compensated though by an increase of e.g. wg, 
as has been shown in Fig. (main text). Only for very short cells smaller than 1.9/im, the model 
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predictions deviate from the experimental data for the parameter choice presented in Fig. 5. In 
this length regime, it predicts a decrease of the mean residence time towards smaller cell lengths. 

G. Analysis of Alternative Models with Regard to Stochastic Switching 

We also analyzed stochastic versions of models presented in (S6) and in (S3) with regard to 
stationary states and stochastic switching in small E. coli cells. The stochastic simulations that we 
performed were one-dimensional and relied on a particle based description together with spatial 
binning (Ax =0.13/xm). The reaction dynamics and particle exchange between the spatial bins 
due to diffusion was implemented with a Gillespie algorithm (S8). The rates for the attachment, 
detachment and diffusion processes were chosen as indicated in the refs. S3, S6, S7. 

In these models, we did not find switching behavior similar to what we observed in cells. We 
systematically scanned the neighborhood of the parameter space of the models given in (S6, S7) 
and (S3), respectively, including independent increase of the parameters up to a factor of 2 and 
decline up to a factor of 0.5. Particle numbers were varied in the simulations up to 20% from the 
originally given values. We did not chose a broader range, because the particle densities are chosen 
according to experimental findings (S2).We always tested cell lengths between 2 and 6 microns. 
The solutions we found for the models presented in (S3, S6, S7) were either oscillations or noisy 
homogeneous states. 



[1] B. D. Corbin, X. C. Yu, and W. Margolin, EMBO J. 21, 1998 (2002). 

[2] Y. Shih, X. Fu, G. King, T. Le, and L. Rothfield, EMBO J. 21, 3347 (2002). 

[3] M. Loose, E. Fischer-Friedrich, J. Ries, K. Kruse, and P. Schwille, Science 320, 789 (2008). 

[4] F. Tostevin and M. Howard, Phys. Biol. 3, 1 (2006). 

[5] G. Meacci, Min oscillations in Escherichia coli. Physical Aspects: Experiments and Theoretical 

Descriptions (VDM, Saarbriicken, Germany), pp. 67-68, (2009). 

[6] K. C. Huang, Y. Meir, and N. S. Wingreen, Proc. Natl. Acad. Sci. USA 100, 12724 (2003). 

[7] K. C. Huang and N. S. Wingreen, Phys. Biol. 1, 229 (2004). 

[8] D. Gillespie, J. Phys. Chem. 81, 23402361 (1977). 



22 



# of stochasticly switching cells 

# of regularly oscillating cells 



160-1 



120- 



80- 



40- 








1.5 2 2.5 3 3.5 4 

cell length (jim) 

FIG. 7: Histogram of Min switching behavior of cells in a given cell length interval. The red column 
counts stochastically switching cells whereas the blue column represents the regularly switching 
cells. The decision whether a cell switches stochastically or not was taken as follows: residence 
times from a cell where sorted according to the cell length interval which the cell was in at the 
beginning of the residence time. If there were more than 3 residence times associated to a given 
cell length interval, with mean smaller than 80s, than we calculated r- lc = o\ c jf\ c , where f\ c is the 
mean and a lc the standard deviation. If this quantity was greater than one, the cell was decided 
to switch stochastically in this cell length interval. If there were less than 3 residence times in 
the corresponding length interval or if the mean f; c was greater than 80s, then we looked at the 
maximal residence time occurring. If it was greater than 100s, we labeled the cell as stochastically 
switching in the respective length interval. 
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FIG. 8: The evolution of stochasticity of Min dynamics with increasing cell length. Yellow triangles: 
population average of n c , where r lc denotes the ratio of the standard deviation to the mean of 
residence times for an individual cell in a time intervall in which the cell grew less than 0.2/im, 
see text for details. Green triangles: The same data of residence times were first pooled for cells 
in the same length interval. Then, ratio r = <r/(r) of the standard deviation to the mean of these 
residence times were calculated. 
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FIG. 9: Examples of missed switching events in otherwise regularly oscillating cells. The top panel 
shows a kymograph from a real cell. The bottom panel was obtained from a simulation. 
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FIG. 10: Kymograph of stochastic switching of fluorescent labeled MinE in two cells of the E.coli 
strain WM1079. The kymographs cover a time span of 40 minutes. 



# of cells 



5 - 



-40 -20 20 40 % 

Relative deviation of MinD concentration 
from equipartition value in daughter cells 

FIG. 11: Histogram of the deviation from the equipartition value I eq = (I^j + I^)/2 of the inte- 
grated MinD fluorescence lP t in a daughter cell after cell division. Given is the relative deviation 
lint/Ieq i n percent. The red part of the columns indicates the number of daughter cells which 
switched stochastically after division and the blue part of the columns those which oscillated regu- 
larly. Note that, since we assume that no GFP-MinD are lost during division, the overall histogram 
has to be symmetric around zero. 
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FIG. 12: Mean residence time (green dots) and standard deviation (black squares) in a simulated 
cell for increasing attachment rate cue of MinE but at constant cell length. The MinD/MinE ratio is 
kept constant as 8/3. Standard deviation and residence time both fall for increasing concentrations. 
Their ratio also decreases, implicating reduced stochasticity of Min dynamics. 
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